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Abstract 

We extend the Aw-Rascle macroscopic model of car traffic into a two-way multi- 
lane model of pedestrian traffic. Within this model, we propose a technique for the 
handling of the congestion constraint, i.e. the fact that the pedestrian density cannot 
exceed a maximal density corresponding to contact between pedestrians. In a first 
step, we propose a singularly perturbed pressure relation which models the fact that 
the pedestrian velocity is considerably reduced, if not blocked, at congestion. In a 
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second step, we carry over the singular limit into the model and show that abrupt 
transitions between compressible flow (in the uncongested regions) to incompressible 
flow (in congested regions) occur. We also investigate the hyperbolicity of the two- 
way models and show that they can lose their hyperbolicity in some cases. We study 
a diff'usive correction of these models and discuss the characteristic time and length 
scales of the instability. 

1 Introduction 

Crowd modeling and simulation is a challenging problem which has a broad range of 
applications from public safety to entertainment industries through architectural and 
urban design, transportation management, etc. Common and crucial needs for these 
applications are the evaluation and improvement (both quantitatively and qualitatively) of 
existing models, the derivation of new experimentally-based models and the construction 
of hierarchical links between these models at the various scales. 

The goal of this paper is to propose a phenomenological macroscopic model for pedes- 
trian movement in a corridor. A macroscopic model describes the state of the crowd 
through locally averaged quantities such as the pedestrian number density, mean velocity, 
etc. Macroscopic models are opposed to Individual-Based Models (IBM's) which follow 
the location and state of each agent over time. Macroscopic models provide a description 
of the system at scales which are large compared to the individuals scale. Although they 
do not provide the details of the individuals scale, they are computationally more efficient. 
In particular, their computational cost does not depend on the number of agents, but only 
on the refinement level of the spatio-temporal discretization. In addition, by comparisons 
with the experimental data, they give access to large-scale information about the system. 
This information can provide a preliminary gross analysis of the data, which in turn can 
be used for building up more refined IBM's. This procedure requires that the link be- 
tween the microscopic IBM and the macroscopic model has been previously established. 
Therefore, macroscopic models which can be rigorously derived from IBM's are crucial. 

The present work focuses on a one-dimensional model of pedestrian traffic in corridors. 
This setting has several advantages: 

1. It makes the problem essentially one-dimensional and is a preliminary step for the 
development of more complex multi-dimensional problems. The present work will 
consider that pedestrian traffic occurs on discrete lanes. This approximation can 
be viewed as a kind of discretization of the actual two-dimensional dynamics. It 
prepares the terrain for the development and investigation of truly two-dimensional 
models. 

2. We can build up on previous experience in the field of traffic flow models. Our 
approach relies on the Aw-Rascle model of traffic flow |3] , which has been proven an 
excellent model for traffic flow engineering [12]. In the present work, this approach 
will be extended to two-way multi-lane traffic flow of pedestrians. 

3. It is easier to collect well-controlled experimental data in corridors than in open 
space (see for instance [M]). 
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4. The relation of the macroscopic model to a corresponding microscopic IBM is more 
easily established in the one-dimensional setting. In |2], it has been proven that 
the Aw-Rascle model can be derived from a microscopic Follow-the-Leader model 
of car traffic. The proof uses a Lagrangian formulation of the Aw-Rascle model. 
The correspondence between the Lagrangian formulation and the IBM cannot be 
carried over to the two dimensional case, because of the very special structure of 
the Lagrangian model in one- dimension. 

The most widely used models of pedestrian traffic are IBM's. Several families of 
models have been developed. Rule based models [38] have been used in particular for 
the development of games and virtual reality, with several possible levels of description. 
But their aim is more to have a realistic appearance rather than really reproducing a 
realistic behavior. More robust models are needed for example to test and improve the 
geometry of various types of buildings. Physicists have proposed some models inspired 
from the fluid simulation methods. In the so-called 'social force' model [181 1201 121], the 
equations of motion for each pedestrian have the form of Newton's law where the force 
is the sum of several terms each representing the 'social force' under consideration. It 
obviously relies on the analogy existing between the displacement of pedestrians and the 
motion of particles in a gas. It describes quite well dense crowds, but not the individual 
trajectories of a few interacting pedestrians. Other approaches have been developed in 
the framework of cellular automata [9l [161 EI] • these models, the non-local interactions 
between pedestrians are made local through the mediation of a virtual floor fleld. These 
models also are meant to describe the motion of crowds, not of individuals. Besides, a 
systematic study of the isotropy of cellular automata models is still lacking. More recently, 
some geometrical models have been developed. Pedestrians try to predict each others' 
trajectories, and to avoid collisions [171 ESI SO]- The knowledge of other pedestrians' 
trajectories depends on the perception that the pedestrian under consideration has, which 
may vary with time. [3S] takes into account the fact that this knowledge is acquired 
progressively. Another type of perception based on the visual fleld is proposed in [52] • 
These models describe well the individual trajectories of a few interacting pedestrians, 
but it is not obvious yet whether they can handle crowds. 

By contrast to microscopic IBM's, macroscopic crowd models are based on the anal- 
ogy of crowd flow with fluid dynamics. A flrst approach has been proposed in [22]. In 
|19] . a fluid model is derived from a gas- kinetic model through a moment approach and 
phenomenological closures. Recently, a similar approach has been proposed in [1]. In 
[23 } [2 ^ [25], a continuum model is derived through optimal control theory and differential 
games. It leads to a continuity equation coupled with a potential fleld which describes the 
velocity of the pedestrians. Other phenomenological models based on the analogy with 
the Lighthill-Whitham- Richards model of car traffic have been proposed by [H [101 [TT] . 
In [5^ [57] , instead of considering a continuous time evolution described by PDE's, the 
evolution of measures is performed on a discrete time scale. In the present paper, we 
shall consider a continuous time description. Macroscopic models provide a description of 
the system at large spatial scales. They can be heuristically justified for a long corridor 
stretch like a subway corridor, when the spatial inhomogeneities are weak (such as low 
variations of the density or velocity in the direction of the corridor). Of course, they 
cannot be used when the spatial inhomogeneities are at the same scale as for instance 
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the mean-interpedestrian distance in the longitudinal direction. In the case of narrow 
corridors, this mean-interpedestrian distance is larger because there are less pedestrians 
in a cross- sect ion, and the condition of weak spatial inhomogeneities is more stringent. 
From a rigorous standpoint, the derivation of macroscopic models from Individual-Based 
models requires that the number of agents be large, which is obviously questionable in 
most situations in pedestrians and highway traffic. Still there is a large literature devoted 
to macroscopic models which seem to provide adequate models for large scale dynamics. 

We will be specifically interested in two-way multi-lane traffic flow models with a 
particular emphasis on the handling of congestions. These points have been previously 
addressed in [H] for pedestrian counter- flows, [39] for multi-lane traffic and [29l [30] for the 
treatment of congestions. However, to the best of our knowledge, none of these different 
features have been included in the same model at the same time. The most difficult point 
is the treatment of congestions. In the recent approach [2S1 13D] the congestion constraint 
(i.e. the limitation of the density by a maximal density corresponding to contact between 
pedestrians) is enforced by means of convex optimization tools (for IBM's) or techniques 
borrowed from optimal transportation such as Wasserstein metrics (for continuum mod- 
els). However, these abstract methods do not leave much space for parameter fltting to 
data and cannot distinguish between the behavior of pedestrians and say, sheep. Our 
technique relies on the explicit derivation of the dynamics of congestions, in the spirit of 
earlier work for traffic [6], [71 [13] . This procedure was initiated in the seminal work [8] . 

The outline of the paper is as follows. We flrst present the modeling approach for a 
one-way one-lane Aw-Rascle model (IW-AR) of pedestrian flow in corridors in section [2l 
We then successively extend this model into a two-way one-lane Aw-Rascle model (2W- 
AR) in section [3] and to a two-way multi-lane Aw-Rascle model (ML-AR) in section [H In 
each section, we present the corresponding Aw-Rascle model, together with a simplifled 
version of it supposing that the pedestrian desired velocity is constant and uniform. We 
refer these simplifled models as "Constant desired velocity Aw-Rascle" (CAR) models. 
Therefore, we successively have the IW-CAR, 2W-CAR and ML-CAR models as Constant 
desired velocity versions of respectively the IW-AR, 2W-AR and ML-AR models. The 
IW-CAR model can be recast in the form of the celebrated Lighthill-Whitham-Richards 
(LWR) model of traffic. 

Finally, for each of these models, we propose a specific treatment of congestion regions. 
This treatment consists in introducing a singular pressure in the AR model which tends 
to infinity as the density approaches the congestion density (i.e. the density at which the 
agents are in contact to each other). This singularly perturbed pressure relation provides 
a significant reduction of the fiow when the density reaches this maximal density. A small 
parameter e controls the thickness of the transition region. In the limit e — > 0, two phases 
appear: an uncongested phase where the fiow is compressible and a congested phases 
where the fiow is incompressible. The transition between these two phases is abrupt, by 
contrast to the case where e stay finite, where this transition is smooth. The location 
of the transition interface is not given a priori and is part of the unknowns of the limit 
problem. 

Table [1] below provides a summary of the various proposed models and their relations. 
One interesting characteristics of two-way models as compared to one-way models 
is that they may lose their hyperbolicity in situations close to the congestion regime. 
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Table 1: Table of the various models, with some of their characteristics and the sections in 
which they are introduced. The meaning of the acronyms is as follows: AR='Aw-Rascle model', 
CAR='Aw-Rascle model with Constant Desired Velocity'. The left column (basic model) refers 
the general formulation of the model and the middle and right columns, the modified models 
taking into account the congestion phenomena. The middle column corresponds to a smooth 
transition from uncongested state to congestion while the right column corresponds to an abrupt 
phase transition. 



Although, this loss of hyperbolicity can be seen as detrimental to the model, the resulting 
instability may explain the appearance of crowd turbulence at high densities. We note 
that a loss of (strict) hyperbolicity has already been found in a multi-velocity one-way 
model [5]. In order to gain insight into this instability, in section [5l we analyze the 
diffusive perturbation of the two-way Aw-Rascle model with constant desired velocity, 
and exhibit the typical time scale and growth rate of the so-generated structures. These 
observables can be used to assess the model and calibrate it against empirical data. In 
order to illustrate these considerations, we show numerical simulations that confirm the 
appearance of these large-scale structures which consist of two counter-diffusing crowds. 
In these simulations, which are presented for illustrative purposes only, in order to explore 
what kind of structures the lack of hyperbolicity of the model leads to, we assume a smooth 
pressure-density relation. Thanks to this assumption, we omit to treat the congestion 
constraint, which is a difficult stiff problem, for which special methods have to be designed 
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(see e.g. [151 [H] ^oi the case of Euler system of gas dynamics and [6l[T3] for the AR model). 



2 One-way one-lane traffic model 

2.1 An Aw-Rascle model for one-lane one-way pedestrian traffic 

In this section, we construct a one-lane one-way continuum model of pedestrian traffic in 
corridors. In this model, we will pay a particular attention to the occurrence of conges- 
tions. We encode the congestion effect into a constraint of maximal total density. This 
work is inspired by similar approaches for vehicular traffic, which have been developed in 



For that purpose, the building block is a one-lane, one-way Aw-Rascle (IW-AR) model 
which has been proposed for vehicular traffic flow [3]. This model belongs to the class of 
second-order models in the sense that it considers that both the density and the velocity 
are dynamical variables which are subject to time-differential equations. By contrast, flrst 
order models use the density as the only dynamical variable and prescribe the density flux 
as a local function of the density. The Aw-Rascle model with constant desired velocity 
considered in section 12.21 is an example of a flrst order model. 

Definition 2.1 (IW-AR model) Let p{x,t) G M the density of pedestrians on the lane, 
u e M+ their velocity, w{x,t) G M+ the desired velocity of the pedestrian in the absence 
of obstacles and p{p) the velocity offset between the desired and actual velocities of the 
pedestrian. The 1 W-AR model is written: 



In this model, the offset p{p) is an increasing function of the pedestrian density. By 
analogy with fluid mechanics, this offset will be often referred to as the pressure, but its 
physical dimension is that of a velocity. 

Using the mass conservation equation, we can see that the desired velocity is a La- 
grangian quantity (i.e. is preserved by the flow), in the sense that: 



It is natural, since the desired velocity is a quantity which is attached to the particles and 
should move together with the particles at the flow velocity. 

This model has been studied in great detail in [3] and proven to derive from a foUow- 
the- leader model of car traffic in [2]. Of particular interest is the fact that this model is 
hyperbolic, with two Riemann invariants. The flrst one is obviously the desired velocity w 
as (jl]) testifles. The second one is less obvious but is nothing but the actual flow velocity 
u. Indeed, from (jl]) and using (fl]), we get: 



ig Ellis]. 



dtp + d:,{pu) = 0, 
dt{pw) + dxipwu) = 
w = u + p{p)- 



(1) 
(2) 
(3) 



dfW + ud^w = 0. 



(4) 



dfU + udxU 



-{dtp + ud^p) 
-p{p){dtP + ud^p) 
p'{p)pd^u, 
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and therefore 



dtu +{u- p\p)p)dxU = 0. (5) 

Therefore, information about the fluid velocity propagates with a velocity 

Cu = u-p{p)p. (6) 

In the reference frame of the fluid, this gives raise to waves moving upstream the flow 
with a speed equal to —p'{p)p. 

Remark 1 We can also consider the evolution of pu instead of that ofu. We obtain from 
0) and using 

dt {pu) + dx {puu) = -{dt {pp) + i9x {ppu) ) 
= -p{dtp + ud^p) 
dp 

where we have introduced the material derivative d/dt = dt + udx- This form is motivated 
by the observation that drivers do not react to local gradients of the vehicle density 
but rather to their material derivative in the frame of the driver. This modification to 
standard gas dynamics like models of traffic was crucial in obtaining a cure to the various 
deficiencies of second order models as observed by Daganzo fW^ . Eq. ^ can also be put 
in the form 

dt{pu) + dx{puw) = -dt{pp) 

= -n\p)dtp 

= n'{p)d,{pu), (8) 



with 



7r(p) = PPip), 7r'(p) = pp'ip)+pip). (9) 



We will consider the IW-AR model as a building block for the pedestrian model. In 
order to make the connection with a microscopic view of pedestrian flow, we consider a 
subcase of this model in the section below. 



2.2 Constant desired velocity 

This one-way Constant Desired Velocity Aw-Rascle (IW-CAR) model assumes that the 
pedestrians can have only two velocities: either a fixed uniform velocity V which is the 
same for all pedestrians and does not vary with time ; or zero, indicating that they 
are immobile. In other words, if because of the high density of obstacles in front, the 
pedestrians cannot proceed further with the velocity V, they have to stop. 
In this case, 

w = V 
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UiP) 




Figure 1: Density flux f{p) = p{V — p{p)) as a function of p in the IW-CAR model. 

is a fixed value and therefore, the actual flow velocity 

u = V-p{p) (10) 

is a local function of p. This leads to a first-order model where the flux velocity is given 
as a local prescription of the density. 

Definition 2.2 (IW-CAR model) Let p{x,t) G M the density of pedestrians on the 
lane, V G M+ the (constant) desired velocity of pedestrians and p{p) the pressure. The 
IW-CAR model is written: 

dtp + dMv-pip))) = o. (11) 

We denote by /(p) = piV — p{p)) the mass flux. The quantity p{p) being an increasing 
function of p, f{p) has a concave shape (and is actually concave if pp{p) is convex), which is 
consistent with classical first-order traffic models such as the Lighthill-Whitham- Richards 
(LWR) model [28]. Figured] provides a graphical view of /(p). It is interesting to note that 
the original IW-AR model can be viewed as a LWR model with a driver- dependent flux 
function f{p,w) = p{w—p) where w is the driver dependent parameter, and consequently 
moves with the flow speed. It follows that the LWR is a useful lab to test concepts 
ultimately applying to the IW-AR model. However, some of the features of the LWR 
model are too simple (such as the conservation of the maxima and minima of p) and a 
realistic description of the dynamics requires more complex models such as the IW-AR 
model. 

It is also instructive to write the IW-CAR model as a second order model, like the 
IW-AR model. Indeed, using ^ and ffTU]) . we can write ffTTl) as: 

dtp + d,{pu) = Q, (12) 
dt{pu) + d^ipuV) = Tx'{p)d^{pu). (13) 

Conversely, if p and u are solutions of this model, using the fact that V is a. constant 
together with eq. (IT^ to modify the second term of fll3p . and using the r.h.s. of eq. ([H]) 
to modify the r.h.s of (IT^ . we find, : 

dt{p{p + u-V)) = Q. 

Therefore, if ffTOl) is satisfied initially, it is satisfied at all times and we recover f|TT]) . 



8 



Remark 2 The IW-CAR model in the form [W\] . [W\] has an interesting interpretation 
in terms of microscopic dynamics, when the pedestrians have two velocity states, the 
moving one with velocity V and the steady one, with velocity 0. Indeed, denoting by 
g{x,t) the density of moving pedestrians and by s{x,t) that of steady pedestrians, we have 

p = g + s. 

Because the moving pedestrians move with velocity V , we can write the pedestrian flux pu 
as 

pu = Vg. (14) 

Since by (EP, pu = p{V — p{p)), we deduce that 

Pip)^ pip) 
9 = pil-^), ^ = P-^- 

Not surprisingly, the offset velocity scaled by the particle velocity is nothing but the pro- 
portion of steady particles and it is completely determined by the total density p. 
We deduce from [l4\ ) that system [W\) . 17^) can be rewritten in the form: 

dt{g + s) + d.,{yg) = Q, (15) 
dt{Vg) + d^{V^g) = Vn\p)d^g, (16) 

Dividing ( fi^) by V and subtracting to / I73]) . we find: 

dtg + d^{Vg) = 7r'{p)d^g, 
dts = -7r'{p)dr,g. 

Thus, the term 7i'{p)dxg represents the algebraic transfer rate from immobile to moving 
particles, while —'n:'{p)dxg represents the opposite transfer. Therefore, this model assumes 
that the pedestrians decide to stop or become mobile again, based not only on local ob- 
servation of the surrounding density, but on the observation of their gradients. More 
precisely, keeping in mind that n'{p) and V have the same sign, the transfer rate from 
the immobile to moving state is positive if the moving particle density increases in the 
downstream direction, indicating a lower congestion. Symmetrically, the transfer rate 
from the moving to immobile state increases if the moving particle density decreases in 
the downstream direction, indicating an increase of congestion. These evaluations of the 
variation of the moving particle density derivative are weighted by increasing functions of 
the density, meaning that the reactions of the pedestrians to their environment are faster 
if the density is large. 

We now turn to the introduction of the density constraint in the IW-AR or IW-CAR 
models. 
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2.3 Introduction of the maximal density constraint in the IW- 
AR model 

The maximal density constraint (also referred to below as the congestion constraint) is 
implemented in the expression of the velocity offset or pressure p. Two ways to achieve 
this goal are proposed. 

In the first one, p is a smooth function of the particle density which blows-up at the 
approach of the maximal allowed density p* . 

In the second one, congestion results in an incompressibility constraint which produces 
non-local effects with infinite speed of propagation of information. In congested regions, 
the pressure is no longer a function of the density but becomes implicitly determined by 
the incompressibility constraint. The transition from uncongested to congested regions is 
abrupt and appears as a kind of phase transition. This second approach can be realized 
as an asymptotic limit of the first approach where compression waves (or acoustic waves 
by analogy with gas dynamics) propagate at larger and larger speeds (so-called low Mach- 
number limit). 

Below, we successively discuss these two strategies. Then, we specifically consider the 
introduction of the congestion constraint within the IW-CAR model. 

2.3.1 Congestion model with smooth transitions between uncongested and 
congested regions 

To implement the congestion constraint, we will highly rely on previous work [6l [TJ [13], 
where this constraint has been implemented in the IW-AR model. We take a convex 
function p{p) such that j9(0) = 0, p'(0) > and p(p) — oo as p — )■ p*. More explicitly, we 
can choose for instance for the pressure: 

p{p)=f{p) = P{p) + Q'{p), (17) 

P(p) = Mp'", m > 1, (18) 

Q'{p)= 7>1. (19) 

p p* 

P(p) is the background pressure of the pedestrians in the absence of congestion (and 
is taken in the form of an isentropic gas dynamics equation of state). is a correction 
which turns on when the density is close to congestion (i.e. e <C 1 is a small quantity), and 
modifies the background pressure to have it match the congestion condition p(p) — )■ oo as 
p^p*. 

Indeed, as long as p — p^, is not too small, the denominator in ( fT9ll is finite and Q^(p) is 
of order e. Thus the pressure p is dominated by the P term. However, a crossover occurs 
when 

1 1 

. ) ~ 

p p" 

i.e. when 

p* - p ~ pp*e^/^. (20) 
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Figure 2: Schematic representation of Q^, -P, and = P + Q^, as a function of p. 

Thus in a density range near p* which scales as e^l^ ^ the correction Q^{p) becomes of 
order unity. 

This is represented schematically on Figure O Note that the precise shape of the term 



narrow region close to congestion. The chosen expression ensures that Q'^{p = 0) = 0, 
and that it becomes significant in the vicinity of p* only. Note also that is an increasing 
function of p, in order to keep the problem hyperbolic. 

The pressure singularity at p = p* ensures that the congestion density p* cannot be 
exceeded. Indeed, let us consider a closed system (e.g. the system is posed on an interval 
[a, b] with periodic boundary conditions) for simplicity. Let uq and wq be the initial 
conditions and suppose that they satisfy 



for some constants u^, um, w^, wm- Then, [3] notices that, at any time, u and w satisfy 
the same estimates: 

< Mm < u{x, t) < Um, < Wm < w{x, t) < Wm, V(x, t) e [a, b] x ]R+. (21) 

In other words, this estimate defines an invariant region of the system. It follows from 
the fact that, u and w being the two Riemann invariants, they are transported by the 
characteristic fields (see eqs. ([5]), (j4])) and therefore, satisfy the maximum principle. From 
( ETj) . we deduce that w — u = p{p) < wm — Um, and we also have that p{p) > at all 
times. Let be the inverse function of p. Since p maps [0,p*) increasingly to M+, then 
p~^ maps increasingly M+ onto [0,p*), from which the estimate p < p~^{wm — ^m) < p* 
follows. This indeed shows that the constraint p < p* is satisfied at all times. From the 
estimate f l2T|) . we also see that u cannot become negative, so that the estimate w > p{p) 
is also satisfied at all times. 

With this e-dependent pressure, the IW-AR model becomes a perturbation problem. 




is not important, as it does not contribute to the pressure law, except in a 



< Mm < Mo < Mm, < M;m < M7o < Wm, 



11 



written as follows: 



dtip'w') + d^ip'w'u') = 0, 



(22) 
(23) 
(24) 



The next section investigates the formal e — )■ limit. 

2.3.2 Congestion model with abrupt transitions between uncongested and 
congested regions 

In the limit e — )■ 0, the uncongested motion remains unperturbed until the density hits 
the exact value p*. Once this happens, congestion suddenly turns on and modifies the 
dynamics abruptly. In the uncongested regions, the flow is compressible ; it becomes 
incompressible at the congestion density p*. Therefore, in the limit e — >■ 0, the abrupt 
transition from uncongested motion (when p < p*) to congested motion (when p = p*) 
corresponds to the crossing of a phase transition between a compressible to an incom- 
pressible flow regime. 

In the limit e — )■ 0, the arguments of [6], [71 [13] can be easily adapted. Suppose that 
— p < p*. In this case, Q^(p^) — )■ and we recover an IW-AR model associated to the 
pressure -P(p): 



If on the other hand, p^ — j- p*, then Q^{p'^) — )■ Q with < Q < wm- Therefore, the 
total pressure is such that p^{p'^) — ?■ p with P{p*) < p- In this case, the model becomes 
incompressible: 



Note that in this congested region, the density does not vary (it is equal to p*) and cannot 
determine the pressure anymore. Indeed, the functional relation between the density and 
the pressure is broken and p may be varying with x even though p does not. The spatial 
variations of p compensate exactly (through f[30|) ) the variations of U7°, in such a way 
that all the pedestrians, whatever their desired velocity is, move at the same speed in the 
congestion region. 

This can also be seen when taking the limit e — )■ in ([S]). Indeed, if p^ — ?■ p* with 
p^{p^){= w^—u^) staying finite, then p^—p* = 0{e^^^) (see ([20!) ) and dp^ /dp ~ e~^^^ — t- oo. 
Therefore, in the congested regime, the derivative of the pressure with respect to the 
density becomes infinite. Inserting this in ([5]) shows that d^u'^ — )■ 0. This ensures that 
all the pedestrians move at the same speed. Simultaneously, this blocks any further 



V + 5x(pV) = o, 
dt{p'w') + d,{p'w\') = Q, 



(25) 
(26) 
(27) 



= 0, 

dtw^ + M°(9^.w° = 0, 

= M° + p, with -P(p*) < p- 



(28) 
(29) 
(30) 
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increase of the density, which cannot become larger than p* . Indeed, the mass conservation 
equation ([1]) tells us that 

= dtp"" + u^^^p^ = -p^d^u^, 

at 

and consequently, if dxU'^ — ?■ 0, any further increase of the density is impeded. 

In the general case, we expect that the two limit regimes coexist. The congested 
region may appear anywhere in the flow, depending on the initial conditions. Congestion 
regions must be connected to uncongested regions by interface conditions. Across these 
interfaces, p and pw, which are conserved quantities obey the Rankine-Hugoniot relations. 
The quantity w, which is thought of as the (locally averaged) pedestrians' desired velocity 
is modified across the interfaces through these relations. However, the bounds fl2T]) are 
preserved (see [3]). 

Connecting congested and uncongested regions is a delicate problem which has been 
investigated in [0] by a careful inspection of Riemann problem solutions. Specifically, [B] 
treats the special case M = in (fT7|) - ffT^ . The present choice of the pressure ffT71) - ffT91) 
is slightly different: in the limit e — > 0, it produces a non-zero pressure in the uncon- 
gested region, while |6] considers that uncongested regions are pressureless in this limit. 
Pressureless gas dynamics develops some unpleasant features (such as the occurrence of 
vacuum, weak instabilities, and so on). Keeping a non-zero pressure in the uncongested 
region in the limit £ — )■ allows to bypass some of these problems and represents an 
improvement over [6]. Of course, the precise choice of m and M must be fitted against 
experimental data. 

We do not attempt to derive interface conditions between uncongested and congested 
regions for the present choice of the pressure. Indeed, the perturbation problem fl2^ - fl2^ . 
even with a small value of e is easier to treat numerically than the connection problem 
between the two models (!25|) - (!27|) and (l28l) - (!30l) . Therefore, we will not regard the limit 
model as a numerically effective one, but rather, as a theoretical limit which provides some 
useful insight. Still, the numerical treatment of the perturbation problem requires some 
care. Of particular importance is the development of Asymptotic-Preserving schemes, i.e. 
of schemes that are able to capture the correct asymptotic limit when e — )■ 0. This is not 
an easy problem because of the blow up of the pressure near p*. Indeed, due to the blow 
up of the characteristic speed in ([5]), the CFL stability condition of a classical explicit 
shock-capturing method leads to a time-step constraint of the type At = 0{e^/^) — )■ 
as £ — 0. For this reason, classical explicit shock- capturing methods cannot be used to 
explore the congestion constraint when e — )■ and Asymptotic-Preserving schemes are 
needed. 

Another reason for considering the perturbation problem (!22|) -( l24|) instead of the limit 
model is that the congestion may appear gradually rather than like an abrupt phase 
transition from compressible to incompressible motion. In particular, for large pedestrian 
concentrations, some erratic motions occur (this is referred to as crowd turbulence) and 
might be modeled by a suitable (may be different) choice of the perturbation pressure Q^. 
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2.3.3 Introduction of the congestion constraint in the constant desired ve- 
locity IW-CAR model 

(i) Congestion model with smooth transitions. The smooth pressure relations f|T71) - f|T9l) 
can be used for the IW-CAR model. Because p now satisfies a convection equation: 

dtp + {v-{ppyip))d,p = o, 

the initial bounds are preserved. Indeed, suppose that 

< Pm < PO < PM < P*, 

for some constants pm, pu, then, at any time, p satisfy the same estimates. 

< Pm < p{x, t) < pm< p*, V(a;, t) e [a, b] x M+. 

In this way, the constraint < p < p* is always satisfied. However, the fact that the 
bounds on the density are preserved by the dynamics can be viewed as unrealistic. In 
real pedestrian traffic, strips of congested and uncongested traffic spontaneously emerge 
from rather space homogeneous initial conditions. The generation of new maximal and 
minimal bounds is an important feature of real traffic systems which is not well taken 
into account in the IW-CAR model and more generally, in LWR models. 

(a) Congestion model with abrupt transitions. If the limit e — )■ is considered, and if the 
upper bound pm = Pm depends on e and is such that p^ — >■ p*, then, some congestion 
regions can occur. The limit model in the uncongested region does not change, and is 
given by the single conservation relation ( |TT|) with the pressure p(p°) = -P(p°). In the 
congested region, we have p^ = p*, which implies dxU^ = 0. In terms of the moving and 
steady pedestrian densities, the congested regime means that 

d,g' = 0, s' = p*- g', 

i.e. both the steady and moving pedestrian densities are uniform in the congested region. 

3 Two-way one-lane traffic model 

3.1 An Aw-Rascle model for two-way one-lane pedestrian traffic 

The extension of the IW-AR model to 2-way traffic, denoted by 2W-AR model, may seem 
rather easy, the 2-way traffic is written as a system of two 1-way models. However, we 
will see that the mathematical properties of the 2-way models are rather different from 
their one-way counterpart. 

Definition 3.1 (2W-AR model) Let p± the density of pedestrians, u± their velocity, 
w± their desired velocity andp the pressure, with an index + for the right-going pedestrians 
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and — for the left-going ones. The 2W-AR model for 2-way traffic is written: 



dtp+ + d^{p+u+) = 0, 
dtp^ + d^ip^u^) = 0, 



dt{p+w+) + d^{p+w+u+) = 0, 
dt{p^w-) + dn,{p-w-u^) = 0, 



w+ = u++p{p+,p^), 
w- = -n_ +p(p_,p+) 



The coupling of the two flows of pedestrians in the 2W-AR model is through the 
prescription of the pressures which are functions of the densities of the two species p+ and 
P-. Our conventions are that the desired velocities w±. and the velocity offsets p{p±, p^) 
are magnitudes, and as such, are positive quantities. The actual velocities m± are signed 
quantities: m+ > for right-going pedestrians and m„ < for left-going pedestrians. 
These conventions explain the different signs in factor of the velocities for fl3ip and fl3ip . 
However, we do not exclude that, in particularly congested conditions, the right-going 
pedestrians may have to go backwards (i.e. to the left) or vice-versa, the left-going 
pedestrians have to go to the right. Therefore, we do not make any a priori assumption 
on the sign of u±. For obvious symmetry reasons, the same pressure function is used for 
the two particles, with reversed arguments. The function p is increasing with respect to 
both arguments since the velocity offset of one of the species increases when the density 
of either species increases. 

Some of the properties of the IW-AR system extend to the 2W-AR one. For instance, 
the desired velocities are Lagrangian variables, as they satisfy: 



Unfortunately, the velocities u+ and u_ do not constitute Riemann invariants any longer 
because of the coupling induced by the dependence of p upon p+ and p_. For this reason 
initial bounds on n+ and u„ are not preserved by the flow, as they were in the case of 
the IW-AR model. Since the velocity offsets p(p+,p_) and p(p_,p+) are not bounded a 
priori, the velocities and m_ can reverse sign when the velocity offsets are large. This is 
expected to reflect the fact that a dense crowd moving in one direction may force isolated 
pedestrians going the other way to move backwards. Of course, such a situation is only 
expected in close to congestion regimes. 

Nonetheless, the evolution of the pedestrian fluxes reflects the same phenomenology 
as in the one-way case, namely that pedestrians react to the Lagrangian derivative of the 
pressure, as shown by the following eqs. (which are the 2-way equivalents of eq. ([7j)): 



dtW+ + U+dri:W+ = 0, 

dtw^ + u^dxW^ = 0. 
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where the material derivatives {d/dt)± = dt + u±dx depend on what type of particles is 
concerned. These equations can also be put in the form (equivalent to ([8]) for the IW-AR 
model): 



p{p+, p„) + p+ dx{p+u+) + 

+ P+ d2p\{p^^p_)d:r{p~U^), 



Pip-,P+)+ P- dip\ 



dx{p-u^ 



P- 92p\^p_^p^)dxip+u+), 



(31) 



(32) 



where we denote by dip and d2P the derivatives of the function p with respect to its first 
and second arguments respectively. This form of the equations will be used below for the 
derivation of the Constant Desired Velocity model. 

The 2W-AR model is not always hyperbolic. Before stating the result, we introduce 
some notations. We define: 



C++ 



dip{p+,p-), 
d2P{p-,P+), 



c+- 

c 



d2p{p+,P-), 
diPip-,p+)- 



(33) 
(34) 



We assume that p is increasing with respect to both arguments, which implies that all 
quantities defined by are non- negative. This assumption simply means that the 

pedestrian speed is reduced if the densities of either categories of pedestrians increase. 
For a given state (p+, w+, p_, w^), the fluid velocities are given by: 



p{p+,p^), 

We also define the following velocities 

Cu+ = U+ ~ P+C++, 



u- = -w-+p{p-,p+). 



u_ + p^c. 



These are the characteristic speeds of the IW-AR system. Specifically, is the wave 
at which information about velocity would propagate in a system of right-going pedestri- 
ans without coupling with the left-going ones. A similar explanation holds symmetrically 
for c„_ . 

We now have the following theorem, the proof of which is elementary and left to the 
reader. 

Theorem 3.2 The 2W-AR system is hyperbolic about the state {p+,w+, p-,W-) if and 
only if the following condition holds true: 



A : = 



(35) 



The quantities u±. are two characteristics velocities of the system. If condition is 
satisfied, the two other characteristic velocities are 



A- 



+ c„_ ± V A 



(36) 
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Non-hyperbolicity occurs when the two characteristic velocities c^^ and c„_ of the uncou- 
pled systems are close to each other. In this case, the first term of A is close to zero and 
does not compensate for the second term, which is negative. These conditions happen in 
particular when both velocities c„_|_ and close to zero, which corresponds to the 

densities where the fluxes p+u+, p-U- are maximal as functions of the densities p+, p_ 
respectively. In particular, in the one-way case with constant speed of figure [H this would 
correspond to the point Pmax- These conditions correspond to the onset of congestion. 
Therefore, instabilities linked to the non-hyperbolic character of the model will develop 
in conditions close to congestion. 

The occurrence of regions of non-hyperbolicity is not entirely surprising. The insta- 
bility of two counter-propagating flows is a common phenomenon in fluid mechanics. In 
plasma physics, the instability of two counter-propagating streams of charged particles 
is well known under the two-stream instability. The situation here is extremely similar, 
in spite of the different nature of the interactions (which are mediated by the long-range 
Coulomb force in the plasma case). 

The occurrence of a non-hyperbolic region is often viewed as detrimental, because in 
this region, the model is unstable. On the other hand, self-organization phenomena like 
lane formation or the onset of crowd turbulence cannot be described by an everywhere 
stable model. For instance, morphogenesis is explained by the occurrence of the Turing 
instability in systems of diffusion equations. Here, diffusion is not taken into account and 
the instability originates from a different phenomenon. However, in practice, some small 
but non-zero diffusion always exists. This diffusion damps the small scale structures but 
keeps the large scale structures growing. The typical size of the observed structures can 
be linked to the threshold wave-number below which instability occurs. 

Numerical simulations to be presented in a forthcoming work will allow us to determine 
whether the phenomena which are observed in dense crowds may be explained by this 
type of instability. In section [5l a stability analysis of a diffusive two-way LWR model 
will provide more quantitative support to these concepts. 

3.2 The constant desired velocity Aw-Rascle model for two-way 
one-lane pedestrian traffic 

To construct the two-way constant desired velocity Aw-Rascle model (2W-CAR) for two- 
way one-lane pedestrian traffic, we must set 

w+ = W- = V, (37) 

and 

n+ = 1/ -p(p+,p_), u_ = -1/ +p(p_,p+). (38) 
This leads to the following model: 

Definition 3.3 (2W-CAR model) Let p+ and p_ the densities of pedestrians moving 
to the right and to the left respectively, V the (constant) desired velocity of pedestrians 
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and p the pressure term. The 2W-CAR model is written: 

atp+ + a.(p+(y-p(p+,p_))) = o 

dtp- - d^{p_{V - p{p_, p+))) 



0. 



(39) 
(40) 



These are two first-order models coupled by a velocity offset which depends on the 
two densities. 

We can find the same interpretation of this model in terms of moving and steady 
particles as in the one-way model case. Using fl3T|) . fl32l) and fl371) . we can write: 



dtp+ + d^{p+u+) = 0, 
dtp- + d^{p-u-) = 0, 
dt{p+u+) + d^{p+u+V) = p{p+, p_) 

dt{p-u-) - d^{p-U-V) = - p(p_,p^ 



P+ dip\ 



(p+,p-) 



dx{p+u+) + 



-p+ d2p\(^p^^p_)d^{p-u. 



+ dip\ 
-p- d2p\ 



dxip-u^ 
dx{p+u+) 



(41) 
(42) 

(43) 
(44) 



Conversely, if p+, u+, p_, m_ are solutions of this model, using the same method as in the 
one-way case, we easily find that: 

dt{p±ip±Ui-V)) = 0. 

Therefore, if (155]) is satisfied initially, it is satisfied at all times and we recover fHUjl . 

Now, we denote by g±{x, t) the density of the moving particles and by s±{x, t) that of 
the steady particles with a -|- (respectively a — ) indicating the right-going (respectively 
left-going) pedestrians. Although steady, the pedestrians have a desired motion either to 
the right or to the left, and we need to keep track of these intended directions of motions. 
We have 



We deduce that 



P±= g± + s± and p±u± = ±Vg±. 

p{p-,p+) 



P+ 



Pip+,P- 
V 



P- 



V 



Therefore, the offset velocities p{p+,p-) and p{p-,p+) scaled by the particle velocity V 
represent the proportions of the steady particles s'^ / p+ and / p- respectively. Now, we 
can rewrite (HT]) - (H1|) as follows: 



dt{g+ + s+) + d,{Vg+) = 0, 
dt{g- + s-)-d,{Vg-) = 0, 

dt{Vg+) + d.AV'g+)= [p(p+,p-) 
dt{Vg-)-d,{V'g-) = -\p{p-,p, 



P+ dip\ 



dx{Vg, 



-P+ 92p\^p^^p_)d^{Vg. 



P- dip\ 



{p-,p+) 



(P-.P+) 



dxiVg-) + 

dAVg+). 



By simple linear combinations, this system is equivalent to 



dtg+ + d:,{Vg+) = p{p+,p-) + p+ d^g+ 

dtg- - d^{Vg-) = - p{p-,p+) + p- d^g. + 

+P- d2p\^p_^p^)d:,g+, 
dts+ = - p{p+,P^) + P+ dxg+ + P+ d2p\io^_n) d^g 



dts^ = [p|(p_,p+) +P-9ip|(p_,p+)] d^g^ - p^d2p\(p_,p+)d^g+. 

Like in the one-way model, we find that the transition rates from the steady to moving 
states or vice- versa depend on the derivatives of the concentrations of moving pedestrians. 
Now, both the left and right going pedestrian total densities appear in the expressions of 
the transitions rates for either species. This is due to the coupling through the pressure 
term, which depends on both densities. 

Like the 2W-AR model, the 2W-CAR model is not always hyperbolic. Using the same 
notations as in the previous section, we have the: 

Theorem 3.4 The 2W-CAR system is hyperbolic about the state if and only if 

condition [3B^) is satisfied. In this case, the two characteristic velocities are given by (EQj- 



This can be seen directly from equations fl39p and fl4Up . once they are put under the form 
dt\ ^+ V f ^d J ^+ 1=0. 



p_ y V P-(^-+ J \ p~ 

We refer to the end of section 13.11 for more comments about this property. 

3.3 Introduction of the congestion constraint in the 2W-AR 
model 

3.3.1 Congestion model with smooth transitions 

It is difficult to make a prescription for the function p. Its expression should be fitted to 
experimental data. Here we propose a form which allows us to investigate the effects of 
congestion. We propose: 

P{p+,P-) =P%P+,P~) = P{p) + Q%P+,P^), with p = p+ + /7_ (45) 
P(p) = Mp™, m > 1, (46) 

Q'(p+,P-) = — T7 7>1- (47) 



g(p^ 



p p 



The rationale for this formula is as follows. First, in uncongested regime, we expect 
that the velocity offsets of the right and left going pedestrians are the same, this common 
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offset being a function of the total particle density. Thus, the uncongested flow pressure 
P given by ( 146|) is a function of p only, and has the same shape as in the one-way case. 
Congestion occurs when the total density p becomes close to p*. Therefore, formula 
( HTIl resembles ( IT^ . except for the prefactor q{p+). With this choice of the pressure, we 
anticipate that the constraint 

p = P+ + p^ < P* 

will be satisfied everywhere in space and time, like in the one-way case. 

The prefactor q{p+) takes into account the fact that the velocity offset for the majority 
particle is smaller than that of the minority particle. Therefore, we prescribe q to be an 
increasing function of p+. For further usage, we note the following formula, which follows 
from eliminating ((1/p) — {l/p*)y between Q^{p+,p_) and Q'^{p^, p+): 

q{p+) Q'{p+, p_) = g(p„) g^(p_, p+). (48) 

It is more convenient to express this formula as 

Q%p+,p-) ^ g(p-) 
Q'{p~,P+) (i{p+y 

remembering that q is an increasing function. This formula states that the velocity offset 
for the right and left-going particles are inversely proportional to the ratios of a (function 
of) the densities. Since q is increasing and taking p_ < p_(_ as an example, we deduce 
that the velocity offset of the right-going particles will be less than that of the left-going 
particles. In other words, the flow of the majority category of pedestrians is less impeded 
than that of the minority one. In order to keep Q^(p±,p^) small whenever p < p*, we 
require that g(p±) = 0(1) when p± < p* ■ Physically relevant expressions of q{p±) can 
be obtained from real experiments. A possible extension, that we will not consider here, 
would be to have different functions g+(p+) and g_(p_). This could model the fact that 
for example, a crowd heading towards a train platform could be more pushy than the one 
going in the opposite direction. 

The 2W-AR model with e-dependent pressure becomes a perturbation problem: 

dtpl + d,{plul) = 0, 
5tpi + 5x.(pinl) = 0, 
dtiplwl) + d,ip%wlu%) = 0, 
dtip'_w'_) + d,{plw'_u'_) = 0, 
wl = uI+p%pI,p'_), 
wt = -ut+p%pt,p%). 

3.3.2 Congestion model with abrupt transitions 

This case corresponds to the formal limit e — > of the previous model. Suppose that 
p^ — p < p*. In this case, <5^(p±,p^) and we recover a 2W-AR model associated to 
the pressure -P(p): 
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5tp° +4(p°X) = 0' 
V- + 5x.(p%°) = o, 

a,(p°O + a,(p>°<) = 0, 

+ = 0, 

=4 + P(pO), ^0 >o, 
u;0 = +P(p°), M°<0. 

If on the other hand, — )■ p*, then Q^(p+,pi) — 0+ and Q'^{p'L,p%) — > Fur- 
thermore, following dUD, Q+ and Q- are related by: 

g(p°)QH. = g(p°)g- (49) 

Therefore, the total pressure is such that p^(p+,pl) p+ and p'^{p^_,p\) p- with 
P(p*) < p± and p+ and p_ related through (149|) (with replaced by p± — P(p*)). 

We stress the fact that Q± and consequently p± are not local function of p° , p° (only 
the ratio Q^/Q^ = qi^p^^J) / q{p%) is a local function of p5j_, p° ). Indeed the value of Q± of 
two different solutions of the model may be different, even if the local values of (p*]., p° ) 
are the same. Therefore, there is no local function of (p*]., p^_) which can match the value 
of Q±. 

Then, in this case, the model becomes: 

dtpl + d,{plul) = Q, 

V- + 5.(p%°-) = 0, 

at(p°O + 9.(p°«) = 0, 
9t(p%°) + 9,.(p%%°) = 0, 

w+=M++p+ with P{p*)<p+, 

= -n° with P{p*) < p_, 

p°+p°=p*, (50) 
q{pl) ip^ - Pip*)) = q{p'^) (P- - P{pl). (51) 

Relations (l50l) and (15T]) furnish the two supplementary relations which allow us to compute 
the two additional quantities p+ and p_. The last relation flSTj) specifies how, at congestion, 
the left and right going pedestrians share the available space. We see that this sharing 
relation depends upon the choice of the function q. Obviously, q is an input of the model 
which must be determined from the experimental measurements. If some flow asymmetry 
must be taken into account (like if one crowd is more pushy than the other one), different 
functions q'+(p+) and g-(p-) can be used. 

This model is a system of first-order differential equations in which the fluxes are 
implicitly determined by the constraint (!50|) . As a consequence of this constraint, the 
total particle flux p^j^vPj^ + p^_vP_ is constant within the congestion region. We note the 
difference between this constrained model and the constrained IW-AR model (see section 
I2.3.2p . In the IW-AR model, there was a single unknown congestion pressure p and 
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a single density constraint p = p*. In the 2W-AR model, there are two congestion 
pressures p+ and which play a similar role in the dynamics of their associated category 
of pedestrians. However, there is still a single density constraint, acting on the total 
density p+ + p_ = p*. The additional condition which allows for the computation of 
the two congestion pressures is provided by the 'space-sharing' constraint fl5T]) . The two 
constraints express very different physical requirements and must be combined in order 
to find the two congestion pressures which, themselves, have a symmetric role. 

3.3.3 Introduction of the congestion constraint in the 2W-CAR model 

(i) Congestion model with smooth transitions. The smooth pressure relations f l45p - fl47p 
can be used for the 2W-CAR model. With this pressure relation, we anticipate that the 
bound p < p* is enforced. 

(a) Congestion model with abrupt transitions. If the limit e — is considered, then, the 
limit model in the uncongested region remains of the same form, i.e. is given by 
( HUj) with the pressure given by p(p°,p°) = -P(p+ + p°). In the congested region, using 
the same arguments as in section [3.3.2[ we find that (p+,p'l) satisfies: 



Again, this model gives rise to a system of first order differential equations in which the 
fluxes are implicitly determined by the constraints ( 152|) . (!53|l . As a consequence of this 
constraint, the total particle flux p+M+ + p^_^u^ (where = V — p±) is constant within 
the congestion region. 

4 Two-way multi-lane traffic model 

4.1 A Two-way multi-lane Aw-Rascle model of pedestrians 

We now consider a multi-lane model to describe the structure of the flow in the cross 
sectional direction to the corridor. The models presented so far considered averaged 
quantities in the cross section of the corridor. However, it is a well observed phenomenon 
that two-way pedestrian flow presents interesting spontaneous lane structures (see e.g. 
[9]), with a preferential side depending on sociological behavior: pedestrians show a pref- 
erence to the right side in western countries, while the preference is to the left in Japan for 
instance. In order to allow for a description of the cross-section of the flow, we discretize 
space in this cross-sectional direction and suppose that pedestrians walk along discrete 
lanes, like cars on a freeway, with lane changing probabilities depending on the state of 
the downwind flow. In this way, we design a model which may, if the parameters are 
suitable chosen, exhibit the spontaneous emergence of a structuration of the flow into 



dtpl + d,{pl{V-p+))=0, 

V_-a.(p°(v^-P_)) = o, 




q{pl) (p+ - P{p*)) = g(P° ) iP- - P{P*))- 



(52) 
(53) 
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lanes. We stress however, that the lanes in our model must be viewed as a mere spatial 
discretization and that spontaneously emerging pedestrian lanes may actually consist of 
several contiguous discrete lanes of our model. 

Let /c G Z be the lane index. So far, we consider an infinite number of lanes. Of course, 
there is a maximal number of K lanes and k G {1, . . . , K}. Extra-conditions due to the 
finiteness of the number of lanes are discarded here for simplicity. For each of the lane, we 
write a 2W-AR model in the form described in section [3l supplemented by lane-changing 
source terms. 

Definition 4.1 (ML-AR model) For any index k G Z, let ^ the density of pedes- 
trians in the k-th lane, Uk^± their velocity, Wk^± their desired velocity and pk a pressure 
term, with an index + for the right-going pedestrians and — for the left-going ones. The 
ML-AR model is given by: 



where Sk,± and Rk,± are source terms coming from the lane-changing transition rates. 

We allow for different pressure relations in the different lanes, to take into account for 
instance that the behavior of the pedestrians may be more aggressive in the fast lanes 
than in the slow ones, or to take into account that circulation along the walls may be 
different than in the middle of the corridor. This point must be assessed by comparisons 
with the experiments. We specify the pressure relation in each lane in the form of fl45p . 
(1471) with parameter values depending on k. 

We denote by 

pk = pk,+ + Pfc,-, 

the total density on the /c-th lane. We assume that the congestion density p* is the same 
for all lanes (this assumption can obviously be relaxed). 

4.2 Interaction terms in the multi-lane model 

We assume that pedestrians prefer to change lane than to reduce their speed, i.e. they 
change lane if they feel that the offset velocity of their lane (i.e. pk{pk,+, pk,-) in the 
case of right-going pedestrians on lane k) increases. If facing such an increase, right- 
going pedestrians change from lane k to lanes ± 1 (not changing their direction of 
motion) with rates X^^f^^^. Similarly, these rates are X'l^^f^^^ for left-going pedestri- 
ans. These rates increase with the value of {d/dt)k,+ {pk{Pk,+ , Pk,-)) for X'l^^f^^^ and with 
{d/dt)k-{pk{Pk,-, Pk,+)) for A^_^^_|_^ to indicate that the lane changing probability is in- 
creased when an increase of the downstream density is detected. We have denoted by 



dtpk,+ + d,j,{pk,+ Uk,+) = Sk,+, 
dtpk,- + d^{pk,- Uk-) = Sk,-, 



(54) 
(55) 
(56) 
(57) 
(58) 
(59) 



dt{pk,+ Wk,+) + dx{pk,+ Wk,+ Uk,+) 
dt{pk,- Wk,-) + dx{pk,- Wk,- Uk,-) 




Wk,+ = Uk,+ +Pk{Pk,+ ,Pk,-), 
Wk,- = -Uk,- +Pk{pk,-,Pk,+)- 
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{d/dt)k^± the material derivatives for particles moving on the /c-th lane in the positive 
or negative direction: {d/dt)k^± = dt + Uk^±dx- Strongly congested lanes do not attract 
new pedestrians. Therefore, A^_j.^^^ is also a decreasing functions of pk+i which vanishes 
at congestion, when p^+i = p* . Similarly, X^^^^j^^ is decreasing with p^+i and vanishes 
at congestion pk+i = p* and X^^j^_i decreases with p^-i and vanishes at congestion 
Pk-i = P*- 

Given these assumptions on the transition rates, the lane-changing source terms for 
the density equations are written: 

Sk,a = Xk+l-^kPk+l,a + K-l^kPk~l,a — iK^k+l + K^k-l)Pk,a, « = ±. (60) 

It is easy to see that this formulation gives: 

fcez 

which implies the balance equation of the total number of particles moving in a given 
direction: 

dtPa + dja = 0, Pa = ^ Pk,a, ja = ^ Pk,aUk,a, a = ±. 

fcez fcez 

Concerning the rates Rk,±, we consider that Wk^± being a Lagrangian quantity, the 
quantities Pk,±Wk,± vary according to the same rates as the densities themselves. Hence, 
we let: 

Rk,a = ^k+l^k Pk+l,a Wk+1 

-iXt^k+l + K^k-l)pk,aWk,a, a = ±. (61) 

The material derivatives of Wk,± satisfy: 



,^ , • dtWk,+ + Uk,+ d^Wk,+ = i.Rk,+ - Wk,+ Sk,+) = 

dt J pk,+ 

,+ Pk+1,+ , \ , \+ Pk-l,+ / X 

= K+i^k — — {wk+i,+-Wk,+) + K^i^k — — (tffc-i,+-wfc,+), 

Pk,+ Pk,+ 

J k~ ' Pk~ 

Pk+l-, \ , \- Pk~l-/ N 

= K+i^k — [Wk+l,--Wk,-) + K^i^k — [Wk-l,--Wk,-). 

Pk- Pk- 

The right-hand sides of these equations are not zero because the arrival of pedestrians from 
different lanes with a different preferred velocity modifies the average preferred velocity. 

4.3 The 'constant desired velocity version' of the two-way multi- 
lane Aw-Rascle model of pedestrians 

To construct the constant desired velocity Aw-Rascle model for two-way multi-lane pedes- 
trian traffic (ML-CAR model), we must set 

Wk,+ = Wk- = V, (62) 
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and 



Uk,+ = V - p{pk,+, Pk,-), Uk, 



V + p{p, 




(63) 



We can check in this case that Sk,± and Rk,± have been defined in a coherent way by (!60l) 
and (16T|) . i.e. that they are such that equations fl54ti55p and fl56H57l) become equivalent. 
The corresponding model is written: 

Definition 4.2 (ML-CAR model) Let pk,± the density of pedestrians in the k-th lane, 
V the constant desired velocity of pedestrian and pk the pressure term. The ML-CAR 
model is given by: 



where Sk,± is given by l[W\) . 

The features of this model are those of the two-way, one-lane CAR model of section 
13. 2[ combined with the features of the source terms Sk,± as outlined in section 14. 2[ 

4.4 Introduction of the congestion constraint in the multi-lane 



4.4.1 Congestion model with smooth transitions 

The prescription for the pressure functions pk are the same as in section 13.3.11 except for 
a possible /c-dependence of the constants, namely: 



With this pressure law, the ML-AR model becomes a perturbation problem. This is 
indicated by equipping all unknowns with an exponent e. This pressure relation can be 
used in the constant desired velocity model of section 14.31 where all particles move with 
the same speed V. 

4.4.2 Congestion model with abrupt transitions 

This case corresponds to the formal limit e — of the previous model. Suppose that 
Pk Pk < P* ■ In this case, Q%{p% p| _) — > and we recover a ML-AR model associated 




ML-AR model 
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to the pressure Pk{pk)'- 



dtipl.wl_) + dM,-<-<-) = 
If on the other hand, p| p*, the model becomes: 

dtpl^ + dM,+ 4,+) = sl+, 
dtpl^ + dM,^ul^) = sl_, 

dtipl^ <d + 9.ipU <+) = 
dtipl,. wi_) + dM,- <- 4,-) = K-, 

w° + = M° + + pk,+ with P{p*) < pk.+, 
wl^^ = -u%^+pk,^ with P{p*) < pk-, 

0,0 * 

Pk,+ + Pk,- = p . 

QkipU) {Pk,+ - Pkip*)) = Qkipl-) (Pk,- - Pkip*)). 

The source terms are unchanged compared to the e > case, and the interpretation of 
the model is the same as in section I3.3.2I Performing the limit e — ?■ in the constant 
desired velocity model of section 14.31 follows a similar procedure and is left to the reader. 



5 Study of the diffusive two-way, one-lane CAR mod- 
els 

In this section, we restrict ourselves to the 2W-CAR model presented in section 13.21 
(i.e. without the introduction of the maximal density constraint), and we investigate the 
stability of a diffusive perturbation of this model. The goal of this section is to show 
that the addition of a small diffusivity stabilizes the large wave-numbers in the region of 
state space where hyperbolicity is lacking. The threshold value of the wave-number below 
which the instability grows can be related to the size of macroscopic structures observed 
in real crowd flows. 

5.1 Theoretical analysis 

We consider the following model which is a slight generalization of the 2W-CAR model: 

dtp+ + da,f{p+, P-) = S dlp+, (64) 
dtp- - dj{p^, p+) = 6 dlp_. (65) 
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Typically, for the 2W-CAR model, f{p+,p-.) = p+{V — p{p+, p-)) but we do not restrict 
ourselves to this simple flux prescription. The assumptions on / are that for flxed p_, 
the function p+ — )■ /(p+,p_) has the bell-shaped curve of flgured], which is characteristic 
of the LWR flux. For flxed p+, the function p_ — )■ /(p+,p_) is just assumed decreasing, 
meaning that the flux of right-going pedestrians is further reduced as the density of left- 
going pedestrians increases. By symmetry, the diffusivities 5 are assumed to be the same 
for the two species of particles. Of course, the diffusivities may depend on the densities 
themselves, in which case they may be different. But we will discard this possibility here. 
We denote by 

C++ = 5i/(p+,p_), c+_ = 92/(p+,p-), 
c_+ = 92/(p_, p+), = 9i/(p_, p+). 

These quantities are related to those deflned in section 13.11 for the 2W-AR model by 

C++ = Cu+, c+_ = -p+c+_, c — = -Cu_, c_+ = -p_c_+. (66) 

With the assumptions on /, we have that c+_ < 0, c_+ < 0, while C++ (resp. c ) 

decreases from positive to negative values when p+ (resp. p_) increases. 

Any state such that (p+, p_) is independent of a; is a stationary solution. We study the 
linearized stability of the system about these uniform steady states. Denoting by (r+, r_) 
its unknowns, the linear system is written: 

dtr+ + c++9x.r+ + c+_(9^r_ = 5 dlr+, (67) 
dtr^ + c^+dxT^ + c d^r^ = Sdlr_. (68) 

We look for solutions which are pure Fourier modes of the form r± = f±expi{C,x — st) 
where f± is the amplitude of the mode, ^ and s are its wave number and frequency. 
Inserting the Fourier Ansatz into fl67j) . f l68|) leads to a homogeneous linear system for 
(f + , f_ ) . This system has non-trivial solutions if and only if the determinant of the linear 
system cancels. This results in a relation between s and C, (the dispersion relation). In 
this analysis, we restrict to ^ G M and are looking for the time stability of the model. We 
denote by A = s/^ the phase velocity of the mode. 

A given mode remains bounded in time, and therefore stable, if and only if the imag- 
inary part of s is non-positive. In the converse situation, the mode is unstable. The 
system is said linearly stable about the uniform state (p+, p_) if and only if all the modes 
are stable for all ^ G M. In the converse situation, the system is unstable, and it is then 
interesting to look at the range of wave numbers G M which generate unstable modes. 
The following result follows easily from simple calculations: 

Proposition 1 (i) Suppose {p^,p_) are such that the following condition: 

A := (5++ + d__f - 4£+_£_+ > 0, (69) 

is satisfied, then the uniform steady state with uniform densities {p^^pJ) is linearly stable 
about {p+,p-). For any given ^ G M, there exist two modes whose phase velocities A+(^) 
are given by 

A±(0 = ^ [£++ - 5- - 2z6^ ± v^] . (70) 
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(ii) Suppose that {p+,p-) are such that l[6^) is not true. Then, the uniform steady state 
with uniform densities {p^,p_) is linearly unstable about {p+,p^). Moreover, we have 



lAI 
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3 a mode such that Im s > (unstable mode). (71) 



The phase velocity is given by 



C++ - c_ 



2i6^±i^/\A\ 



(72) 



We note that if ( 1661) is inserted in ( l69l) . we recover (1351) . Therefore, the addition of diffusion 
does not change the criterion for stabihty or instability. However, in the unstable case, 
all modes are unstable for the diffusion-free model (this would correspond to 5 = in 
(j7T!) ). The addition of a non-zero diffusivity stabilizes the modes corresponding to the 
small scales (large C,)- However, the large scale modes (small ^) remain unstable. We also 
note that, in the stable case, letting the diffusivity go to zero allows us to recover the 
characteristic speed of the diffusion-free model fl5^ . 

For unstable modes, (172|) provides the typical growth rate Ug: it is equal to the positive 
imaginary part of |^|A+ , and given by 



It is maximal for 



1^1 
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Therefore, the typical length scale of the unstable structures is given by the inverse of 
this wave-number: 

because the other modes, having smaller growth rate, will eventually disappear compared 
to the amplitude of the dominant one. These length scale Lg and time scale l/i'g may 
be related to observations and provide a way to assess the model and calibrate it against 
empirical data. 

5.2 Numerical simulations 

In this part, we want to investigate numerically the system f lM|) . flBS]) and in particular we 
are interested in the profile of the solutions whether the system is in a hyperbolic region 
or not. 

With this aim, we first fix a flux function /(p+,p_) defined as: 

f( \ 9{p+ + P-) /7o\ 
/(P+,P-)=P+ \ > (73) 
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where g is a. flux depending on the total density p = + p-. We choose for g a simple 
function increasing on [0, a] and decreasing on [a, 1]: 



X 



< 



X for < X < a 

2a - - 

a aia — x) 

— for a < X < 1 



2 
I 



2 1 



otherwise 



Note that here, in order to keep the simulations simple, we choose a much smoother 
expression for / than the one that was proposed in section 13.31 to enforce the density 
constraint. As a result, the density here can become larger than p* = 1. 



The flux function /(/)+, p 




Figure 3: Left figure: the flux function f{p^,p^) (1731) used in our simulations. Right 
figure: the region of non-hyperbolicity (l69l) of the model, e.g. A < in this region. 



In the following, we take the maximum of g to be at .7, e.g. a = .7. The function / is 
a decreasing function of p_ since g satisfies g'{x) < 1 and / is zero when the total mass 
is greater than 1, e.g. f{p+,p-) = if p+ + p_ > 1. We plot the graph of the function 
/ in figure [3] (left). Then, we numerically compute A to determine the region where the 
system is non- hyperbolic (see figure [3], right). 

To solve numerically the system fl64p . (pHj) . we use a central scheme With this 

aim, we consider a uniform grid in space {xi}i {Ax = Xj+i — Xj) on a fixed interval [0, L] 
along with a fixed time step At. We denote by the approximation of {p^,p-) on the 
cell [xi-i/2,Xi^i/2] (with Xj+i/2 = Xi + Ax/2) at the time nAt. The numerical scheme 
consists of the following algorithm: 



n.+ l 



Xt + A^ - - s , 



Here, -Fi+1/2 denotes the numerical fiux at Xj+1/2 defined as: 



i+1/2 



/2) 



mHl/2) U^Hl/2 
0«+l/2 
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where F is the flux of the system -F(p+, p_) = {f{p+, p_), —f{p-, p+))"^, the vectors 

and ^/^i/2 respectively the left and right value of {p+,p-) at Xj+1/2 computed using a 

MUSCL scheme [27] and aj+i/2 is the maximum eigenvalues fl36l) of the system at and 

Xi+l- 

ai+i/2 = max(|Af I, |A=^i|). 

As initial condition, we use a uniform stationary state {p^,p^) perturbed by stochastic 
noise: 

p+(0,x) = p+ + (Te+(x) , p_(0,x) = p„ + ae_(x), 

with e+(x) and e_(x) two independent white noises and a the standard deviation of the 
noise. We use periodic boundary condition for our simulations. The parameters of our 
simulations are the following: space mesh Aa; = 1, time step At = .2 (CFL= .406), 
diffusion coefficient 6 = A and standard deviation of the noise a = 10^^. We use periodic 
boundary conditions. 

To illustrate our numerical scheme, we use three different initial conditions. First, we 
pick two values for (p+,p_) in the hyperbolic region: 

p+ = .35 , p_ = .3. 

The initial datum is plotted on figure H] (left). As we can see on figure H] (right), the 
solution stabilizes around the stationary state (.35, .3). 

For our second simulation, we take (p+,p_) in a non-hyperbolic region: 

p+ = .5 , p_ = .3. 

The solution does no longer stabilize around the stationary state (.5, .3). On figure [5] 
(left), we observe the apparition of clusters of high density. Each cluster for p+ faces a 
cluster for p_. Moreover, in each cluster, the total mass p+ + p_ is greater or equal to 
1. Therefore the fiux in this region is zero. However, due to the diffusion, the solution is 
not in a stationary state. There is exchange of mass between the clusters. If we run the 
solution for a long time, only one cluster remains (see figure E] (right)). In this cluster, 
we observe that the profile of p_|_ is concave-down whereas the profile of p_ is concave-up. 
Consequently, the diffusivity makes p+ moving backward and p_ moving forward. As a 
result, all the clusters are moving to the left. However, the concavity of the solution is 
puzzling. Numerically, it appears that the concavity of p+ and p_ depends on the total 
mass: the density with higher mass is concave-down and the density with lower mass is 
concave-up. But this property has to be understood analytically. 

For the third simulation, we take an initial datum (p_|_, p_) close to the non- hyperbolic 
region: 

p+ = .4 , p_ = .3. 

Indeed, we can see on figure [3] (right) that the point (0.4, 0.3) almost lies at the border of 
the non-hyperbolic region. The oscillations amplify and clusters of high densities emerge 
(figure El left). However, if we increase the diffusion coefficient, taking 6 = 2 instead of 
6 = .4, then the solution stabilizes around the stationary state (.4, .3) as we observe on 
figure [6], right. Therefore, a large enough diffusion prevents cluster formation. 
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(p+,p_) at time i = 



1.2 
1 
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0.6 
0.4 
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{p+..P^) at time * = 500 




200 



Figure 4: The initial condition (left figure) and the solution at t 
solution stabilizes around the stationary state (.35, .3). 



500 unit times. The 



(p+,P_) at time ( = 5i)0 (p+.p-) at time ( = i()^ 




X X 



Figure 5: Starting from the initial state (.5, .3), the initial oscillations amplify to create 
clusters (left figure). After a longer time (t = 10"* unit times), only one cluster remains 
(right figure). 



(p+,p_) at time * = son (a = .4) (p+,p-) at time ; = 500 [& = 2) 




50 100 150 200 50 100 150 200 



Figure 6: Starting from the initial state (.4, .3), clusters appears once again (left figure). 
However, if we increase the diffusion coefficient (5 = 2 instead of 5 = .4), the solution is 
stabilizing (right figure). 
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The simulations are provided for ilfustration purpose only, and explore what kind of 
structures the lack of hyperbolicity of the model leads to. The experimental evidence of 
the appearance of clusters is difficult to provide since they must occur (if they occur) 
at very high densities. Experiments in such high density conditions are not possible for 
obvious safety reasons. The observation of real crowds shows that pedestrians can still 
move even at very high densities thanks to the spontaneous organization of the flow into 
lanes. The simple one- dimensional model that is simulated here cannot account for this 
feature. However, we conjecture that cluster formation can be impeded in the multi-lane 
model through the introduction of adequate lane-changing probabilities. 

6 Conclusion 

In this work, we have presented extensions of the Aw-Rascle macroscopic model of traffic 
flow to two-way multi-lane pedestrian traffic, with a particular emphasis on the study of 
the hyperbolicity of the model and the treatment of congestions. 

A first important contribution of the present work is that two-way models may lose 
their hyperbolicity in certain conditions and that this may be linked to the generation 
of large scale structures in crowd flows. Adding diffusion helps stabilize the small scale 
structures and favors the development of large scale structures which may be related to 
observations. We have shown numerical simulations which support this interpretation. 

A second contribution of this work is to provide a methodology to handle the conges- 
tion constraint in pedestrian traffic models. Congestion effects reflect the fact that the 
density cannot exceed a limit density corresponding to contact between pedestrians. We 
have proposed to treat them by a modiflcation of the pressure relation which reduces the 
pedestrian velocities when the density reaches this maximal density. If this modiflcation 
occurs on a very small range of densities, then, the model exhibits abrupt transitions 
between compressible flow (in the uncongested region) and incompressible flow (in the 
congested region). 

Mathematically rigorous proofs that these models respect the upper-bound on the 
total density are left to future works. Their numerical resolution will require the develop- 
ment speciflc techniques such as Asymptotic-Preserving methodologies in order to treat 
the occurrence of congestions. Data learning techniques will then be applied to flt the 
parameters of the model to experimental data. Other possible extensions of this work are 
the development of more complex models such as two-dimensional models, kinetic mod- 
els allowing for a statistical distribution of velocities or crowd turbulence models with 
weak compressibility near congestion. Finally, the derivation of approximate equations 
describing the geometric evolution of the transition interface between the uncongested 
and congested regions would help understanding the dynamics of these interfaces. 

References 

[1] S. Al-nasur & P. Kachroo, A Microscopic-to-Macroscopic Crowd Dynamic model, 
Proceedings of the IEEE ITSC 2006, 2006 IEEE Intelligent Transportation Systems 
Conference Toronto, Canada, September 17-20 (2006). 



32 



[2] A. Aw, A. Klar, A. Materne, M. Rascle, Derivation of continuum traffic flow models 
from microscopic follow-the-leader models, SIAM J. Appl. Math., 63 (2002), 259-278. 
MR1952895 (2003m:35148) 

[3] A. Aw, M. Rascle, Resurrection of second order models of traffic flow, SIAM J. Appl. 
Math., 60 (2000), 916-938. MR1750085 (2001a:35111) 

[4] N. Bellomo, C. Dogbe, On the modelling crowd dynamics: from scaling to second 
order hyperbolic macroscopic models. Math. Models Methods Appl. Sci., 18 (2008), 
1317-1345. MR2438218 (2009g:92074) 

[5] S. Benzoni-Gavage, R. M. Colombo, An n-populations model for traffic flow, Euro- 
pean J. Appl. Math., 14 (2003), 587-612. MR2020123 (2004k:90018) 

[6] F. Berthelin, P. Degond, M. Delitala, M. Rascle, A model for the formation and 
evolution of traffic jams. Arch. Rat. Mech. Anal., 187 (2008), 185-220. MR2366138 
(2008h:90022) 

[7] F. Berthelin, P. Degond, V. Le Blanc, S. Moutari, J. Royer, M. Rascle, A Traffic- 
Flow Model with Constraints for the Modeling of Traffic Jams, Math. Models Methods 
Appl. Sci., 18, Suppl. (2008), 1269-1298. MR2438216 (20101:35233) 

[8] F. Bouchut, Y. Brenier, J. Cortes, J. F. Ripoll, A hierachy of models for two-phase 
flows, J. Nonlinear Sci., 10 (2000), 639-660. MR1799394 (2001j:76109) 

[9] C. Burstedde, K. Klauck, A. Schadschneider, J. Zittarz, Simulation of pedestrian 
dynamics using a 2-dimensional cellular automaton, Physica A, 295 (2001), 507- 
525. arXiv:cond-mat/0102397 

[10] C. Chalons, Numerical approximation of a macroscopic model of pedestrian flows, 
SIAM J. Sci. Comput., 29 (2007), 539-555. MR2306257 (2008a:35188) 

[11] R. M. Colombo, M. D. Rosini, Pedestrian flows and nonclassical shocks. Math. Meth- 
ods Appl. Sci., 28 (2005), 1553-1567. MR2158218 (2006b:90009) 

[12] C. Daganzo, Requiem for second order fluid approximations of traffic flow, Transp. 
Res. B, 29 (1995), 277-286. 

[13] P. Degond, M. Delitala, Modelling and simulation of vehicular traffic jam formation, 
Kinet. Relat. Models, 1 (2008), 279-293. MR2393278 (2009a:90022) 

[14] P. Degond, J. Hua, L. Navoret, Numerical simulations of the Euler system with 
congestion constraint, preprint. arXiv: 1008.4045 

[15] P. Degond, M. Tang, All speed scheme for the low Mach number limit of the Isentropic 
Euler equations, Commun. Comput. Phys., 10 (2011), 1-31. arXiv:0908.1929 

[16] R. Y. Guo, H. J. Huang, A mobile lattice gas model for simulating pedestrian evacu- 
ation, Physica A, 387 (2008), 580-586. 



33 



[17] S. J. Guy, J. Chhugani, C. Kim, N. Satish, M. C. Lin, D. Manocha, P. Dubey, 
Clearpath: Highly parallel collision avoidance for multi- agent simulation, in ACM 
SIGGRAPH/Eurographics Symposium on Computer Animation (SCA), pp. 177-187 
(2009). 

[18] D. Helbing, A mathematical model for the behavior of pedestrians, Behavioral Science, 
36 (1991), 298-310. 

[19] D. Helbing, A fluid dynamic model for the movement of pedestrians, Complex Sys- 
tems, 6 (1992), 391-415. 

[20] D. Helbing, P. Molnar, Social force model for pedestrian dynamics, Physical Review 
E, 51 (1995), 4282-4286. 

[21] D. Helbing, P. Molnar, Self- organization phenomena in pedestrian crowds in: F. 
Schweitzer (ed.) Self-Organization of Complex Structures: From Individual to Col- 
lective Dynamics, pp. 569-577, Gordon and Breach, London (1997). 

[22] L. F. Henderson, On the fluid mechanics of human crowd motion. Transportation 
Research, 8 (1974), 509-515. 

[23] S. Hoogendoorn, P. H. L. Bovy, Simulation of pedestrian flows by optimal control and 
differential games. Optimal Control Appl. Methods, 24 (2003), 153-172. MR1988582 
(2004d:93088) 

[24] R. L. Hughes, A continuum theory for the flow of pedestrians. Transportation Re- 
search B, 36 (2002), 507-535. 

[25] R. L. Hughes, The flow of human crowds, Ann. Rev. Fluid Mech., 35 (2003), 169-182. 

[26] A. Kurganov, E. Tadmor, New high-resolution central schemes for nonlinear con- 
servation laws and convection- diffusion equations, J. Comput. Phys., 160 (2000), 
240-282. MRl 763829 (2001c:65102) 

[27] R. J. LeVeque, "Finite volume methods for hyperbolic problems", Cambridge Uni- 
versity Press (2002). 

[28] M. J. Lighthill, J. B. Whitham, On kinematic waves. I: flow movement in long rivers. 
II: A theory of traffic flow on long crowded roads, Proc. Roy. Soc, A229 (1955), 281- 
345. MR0072605 (17,309e) and MR0072606 (17,310a). 

[29] B. Maury, A. Roudneff-Chupin, F. Santambrogio, A macroscopic crowd motion model 
of gradient flow type. Math. Models Methods Appl. Sci., 20 (2010), 1787-1821. 
MR2735914 

[30] B. Maury, J. Venel, A mathematical framework for a crowd motion model, C. R. 
Acad. Sci. Paris, Ser. I, 346 (2008), 1245-1250. MR2473301 (2009m:91150) 



34 



[31] K. Nishinari, A. Kirchner, A. Namazi, A. Schadschneider, Extended floor field CA 
model for evacuation dynamics, lEICE Transp. Inf. & Syst., E87-D (2004), 726-732. 
arXiv:cond-mat / 0306262 

[32] J. Ondfej, J. Pettre, A-H. Olivier, S. Donikian, A synthetic-vision based steering 
approach for crowd simulation, in SIGGRAPH '10 (2010). 

[33] S. Paris, J. Pettre, S. Donikian, Pedestrian reactive navigation for crowd simulation: 
a predictive approach, Eurographics, 26 (2007), 665-674. 

[34] Pedigree team. Pedestrian flow measurements and analysis in an annular setup, in 
preparation. 

[35] J. Pettre, J. Ondfej, A-H. Olivier, A. Cretual, S. Donikian, Experiment- based mod- 
eling, simulation and validation of interactions between virtual walkers, in SCA '09: 
Proceedings of the 2009 ACM SIGGRAPH/Eurographics Symposium on Computer 
Animation, pp. 189-198 (2009). 

[36] B. Piccoli, A. Tosin, Pedestrian flows in bounded domains with obstacles, Contin. 
Mech. Thermodyn., 21 (2009), 85-107. arXiv:0812.4390 

[37] B. Piccoli, A. Tosin, Time-evolving measures and macroscopic modeling of pedestrian 
flow. Arch. Ration. Mech. Anal., 199 (2010), 707-738. arXiv:0811.3383 

[38] C. W. Reynolds, Steering behaviors for autonomous characters, in Proceedings of 
Game Developers Conference 1999, San Jose, California, pp. 763-782 (1999). 

[39] V. Shvetsov, D. Helbing, Macroscopic dynamics of multi-lane traffic, Phys. Rev. E 
59 (1999), 6328-6339. arXiv:cond-mat/9906430 

[40] J. van den Berg, H. Overmars, Planning time-minimal safe paths amidst unpredictably 
moving obstacles. Int. Journal on Robotics Research 27 (2008), 1274-1294. 

[41] W. G. Weng, S. F. Shena, H. Y. Yuana, W. C. Fana, A behavior-based model for 
pedestrian counter flow, Physica A 375 (2007), 668-678. 

[42] M. Zhang, A non- equilibrium traffic model devoid of gas-like behavior. Transportation 
Res. B 36 (2002), 275-290. 



35 



